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The rainbow truncation of the quark Dyson-Schwinger 
equation is combined with the ladder Bethe-Salpeter equa- 
tion for the meson amplitudes and the dressed quark-photon 
vertex in a self-consistent Poincare-invariant study of the pion 
and kaon electromagnetic form factors in impulse approxima- 
tion. We demonstrate explicitly that the current is conserved 
in this approach and that the obtained results are independent 
of the momentum partitioning in the Bethe-Salpeter ampli- 
tudes. With model gluon parameters previously fixed by the 
condensate, the pion mass and decay constant, and the kaon 
mass, the charge radii and spacelike form factors are found to 
be in good agreement with the experimental data. 

Pacs Numbers: 24.85.+p, 14.40.Aq, 13.40.Gp, ll.10.St 



I. INTRODUCTION 

The light pseudoscalar mesons play an important role 
in understanding low-energy QCD. They are the light- 
est observable hadronic bound states of a quark and an 
anti-quark, and are the Goldstone bosons associated with 
chiral symmetry breaking. Their static properties such 
as the mass and decay constants have been studied exten- 
sively Dynamic properties and scattering observ- 
ables are much less understood theoretically, but there- 
fore not less important to calculate within QCD. In this 
respect, the elastic electromagnetic form factors of the 
pion and kaon are very interesting: the probe is well 
understood, there are accurate data for at low Q 2 
to confront theoretical calculations with, and the charge 
radii r~, r 2 K+1 and r 2 KQ are experimentally known. Cur- 
rently, there are several experiments at JLab to deter- 
mine both the pion and the kaon form factor in the range 
0.5 < Q 2 < 3GeV 2 to better accuracy pLH, which could 
help to discriminate between different model calculations. 

To calculate these form factors, we use an approach 
based on the Dyson-Schwinger equations [DSEs], which 
form an excellent tool to study nonperturbative aspects 
of hadron properties in QCD B. The approach is con- 
sistent with quark and gluon confinement |5|6|] , generates 
dynamical chiral symmetry breaking J7j , and is Poincare 
invariant. It is straightforward to implement the correct 
one-loop renormalization group behavior of QCD and 
obtain agreement with perturbation theory in the pertur- 
bative region. Provided that the relevant Ward identities 
are preserved in the truncation of the DSEs, the corre- 
sponding currents are conserved. Axial current conser- 
vation induces the Goldstone nature of the pions and 



kaons |8j ; electromagnetic current conservation produces 
the correct hadronic charge without fine-tuning. 

We obtain the meson Bethe-Salpeter amplitudes 
[BSAs] and the quark-photon vertex as the solutions 
of respectively the homogeneous and inhomogeneous 
Bethe-Salpeter equation [BSE] in ladder truncation. The 
required dressed quark propagators are obtained from so- 
lutions of the quark DSE in rainbow truncation. Non- 
analytic effects from vector mesons are automatically 
taken into account, because these vector qq bound states 
appear as poles in the quark-photon vertex solution ||. 
We employ a realistic model for the effective quark- 
antiquark coupling that has been shown to reproduce the 
pion and kaon masses and decay constants 0] as well as 
the masses and decay constants for the vector mesons 
p, 4> and K* to within 10% flCfl . The model parameters 
are all fixed in previous work |10[| and constrained only 
by m„, ttik, fir and (qq). The produced pion charge ra- 
dius is within 2% of the experimental value Q. Here, we 
use the same approach, without parameter adjustment, 
to calculate the neutral and charged kaon form factors 
and charge radii; we also extend our previous pion form 
factor calculations || to the spacelike Q 2 -domain antici- 
pated for future JLab |i) data. 

In Sec. II we review the formulation that underlies 
a description of the pion and kaon charge form factors 
within a modeling of QCD through the DSEs. Within 
the impulse approximation, we outline the manner in 
which a ladder-rainbow dynamics for the propagators, 
BSAs and quark-photon vertex preserves the meson elec- 
tromagnetic current. We further indicate the additional 
terms needed for current conservation, if one goes beyond 
rainbow-ladder truncation for the DSEs. In Sec. Ill we 
discuss the details of the model and present our numeri- 
cal results for the form factors. Concluding remarks are 
given in Sec. IV. 



II. PSEUDOSCALAR ELECTROMAGNETIC 
FORM FACTORS 

The 3-point function describing the coupling of a pho- 
ton with momentum Q to a pseudoscalar meson with 
initial and final momenta P± = P ± Q/2 respectively can 
be written as the sum of two terms 



Af(P, Q) = Q a Af a (P, Q) + Q b KT(P, Q) 



(1) 



with Q the electric charge of the (anti-)quark, | for the 
w-quark, and — -g for the d- and s-quarks, and with A aba 
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and A abb describing the coupling of a photon to the quark 
and anti-quark inside the meson respectively. The meson 
form factor is defined as 



K\P 1 Q) = 2P U F{Q 2 ). 



(2) 



and the corresponding charge radius as r 2 = — 6F'(Q 2 ) 
at Q 2 = 0. Analogously, we can define a form factor for 
each of the two terms on the RHS of Eq. ([[]) 



A? b (P ) Q) = 2P u F am (Q*). 



(3) 



Current conservation dictates that each of the form fac- 
tors F all {Q 2 ) and F aha (Q 2 ) are 1 at Q 2 = 0. 



A. Impulse Approximation 

Using dressed quark propagators, bound state BSAs, 
and the dressed gg7-vertex, form factors can be cal- 
culated in impulse approximation. We denote by 
r^(<7, <?'; Q) the quark-photon vertex describing the cou- 
pling of a photon with momentum Q to a quark with 
final and initial momenta q and q' = q — Q respectively 
and flavor a. With this notation, the vertices in Eq. (|l|) 
take the form[] 

Af (P,Q) = 2N C J Tr[S a (q)r ah (q,q + :P_) 
xS b (q + )iT b v (q + , q _;Q)S b (q-)f ah ( q _, q ;P + )] , (4) 

where q = k+\P, q± = k-\P±\Q, P± = P ± ±Q, 

and analogously for A^ ba . The notation J A denotes 
a translationally-invariant regularization of the integral, 
with A the regularization mass-scale, which can be re- 
moved at the end of all calculations by taking the limit 
A — > oo. 5(g) is the dressed quark propagator and 
T ab (q,q'-P) is the meson BSA, with P 2 = -to 2 the on- 
shell meson momentum, and q and q' = q — P the quark 
and anti-quark momenta respectively. 

Both S(q) and T ab (q,q';P) are solutions of their re- 
spective DSEs 

5(p) _1 = Z 2 ij>+ Zi,m{fi) 

rA d 4 q „_ , ,A' 



(2tt)< 



g 2 D^{p-q)- 7li S(q)ri(q,p), (5) 



and 



V'-Uk !,':())= [ A 0^K(p,q;Q) X ab (q,q';Q), (6) 



1 We use Euclidean metric {7^,7^} = 25^, 7^ = 7 M and 
a ■ b = J2t=i aibl - 



where Du U (k) is the renormalized dressed-gluon propa- 
gator, T[,(q,p) is the renormalized dressed quark-gluon 
vertex, K is the renormalized qq scattering kernel that 
is irreducible with respect to a pair of qq lines, and 
X al (q,q';Q) = S a (q)T ab (q,q';Q)S b (q') is the BS wave 
function. 

The solution of Eq. (||) is renormalized according to 
5(p) _1 = ip" + m{p) at a sufficiently large spacelike /j, 2 , 
with to(u) the renormalized quark mass at the scale //. 
In Eq. (||), 5, and m(/i) depend on the quark flavor, 
although we have not indicated this explicitly. The renor- 
malization constants Z2 and Z4 depend on the renormal- 
ization point and the regularization mass-scale, but not 
on flavor: in our analysis we employ a flavor-independent 
rcnormalization scheme. 

The meson BSAs T ab (q,q';P) are normalized accord- 
ing to the canonical normalization condition 



Pu 



Nr. 



d 



d 4 q 



dP^J (2tt) 4 



Tr 



r ab (q',q;Q) 



x S a (q + V P) T ab (q, q'; Q) S b {q + (r, - 1)P) 
A d 4 k 



(2nY 



rTr 



X ab (k',k;Q)K(k,q;P) X ab (q,q';Q) 



(7) 



at the mass shell P 2 = Q 2 = — to 2 , with q = q + rjQ, 
q' = q+ {r) — 1)Q, and similarly for k and k! . We use the 
conventions where f n — 92 MeV, and rj describes the mo- 
mentum partitioning between the quark and anti-quark. 
Note that physical observables should be independent of 
this parameter. 

For pseudoscalar bound states the BSA is commonly 
decomposed into 

T(fc + V P, k+( V - 1)P; P) = 
= 75 [iE(k 2 ; k-P;rj)+ PF(k 2 ; k ■ P\rj) 

+ ft G(k 2 ; k-P; V )+ tr M „ k„P v H{k 2 ; k ■ P; v )] , (8) 

with the invariant amplitudes E, F, G and H being 
Lorentz scalar functions of k 2 and k ■ P = kP cos 9. Sub- 
sequently, each invariant amplitude can be expanded in 
k ■ P based on Chebyshev polynomials 

00 

/ (fc 2 , k-P;P 2 )=J2 Ui(co S e)(kPyfi(k 2 ;P 2 ) . (9) 

i=0 

For charge-parity eigenstates, such as the pion, the in- 
variant amplitudes E, F, G, and H have a well-defined 
charge-parity if one chooses rj = 4. Therefore, these am- 
plitudes are either entirely even {E, F, and H) or odd 
(G) in k-P, and one needs only the even (or odd) Cheby- 
shev moments to completely describe these amplitudes, 
which makes this a convenient decomposition. 
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B. The quark-photon vertex 

The quark-photon vertex is the solution of the renor- 
malizcd inhomogeneous BSE 



/A 14 
(2^ K(p ' q > 



Q) 



(2nY 

xS a (q + )r*( q+ ,q^Q)S a (q^, (10) 

with p± — p ± and q± = q ± \Q-, and with the 
same kernel K as the homogeneous BSE for meson bound 
states. Because of gauge invariance, it satisfies the Ward- 
Takahashi identity [WTI] 



(11) 



Solutions of the homogeneous version of Eq. ([Toj ) at 
discrete timelike momenta Q 2 define vector meson bound 
states with masses my = — Q 2 . It follows that Q) 
has poles at those locations, and behaves like 



r>+, P -;Q) 



Tf v ( P+ ,p-;Q)fvm v 



Q 2 



(12) 



in the vicinity of these bound states, where v is the 
aa vector meson BSA, and fy the electroweak decay con- 
stant @. 

For the photon coupled to u- and d-quarks, this re- 
sults in a p-meson pole at Q 2 = —0.6 GeV 2 . For the pho- 
ton coupled to s-quarks, the first pole is located around 
Q 2 = —1.0 GeV 2 at the </>-mass. At the level of the lad- 
der approximation, which is commonly used in practical 
calculations, there is no width generated for the vector 
meson, and the vertex has real poles. One would have to 
incorporate the open tttt channel in the ladder BSE ker- 
nel to produce a vector meson width; for the vertex, this 
would generate an imaginary part beyond the threshold 
for pion production, Q 2 < —4 m 2 in the timelike region. 

The full vertex T° can be decomposed into 4 longi- 
tudinal components and 8 transverse components. The 
longitudinal components do not contribute to the form 
factors. In Ref. || it was shown that only 5 of the 8 
transverse components are important for the pion form 
factor, in the momentum range —0.3 < Q 2 < 1.0 GeV 2 
the remaining 3 components contribute less than 1%. We 
expect that this will also be the case for the kaon form 
factor, and use the Dirac amplitudes Ti to T 5 of Ref. ||] 
only. 



C. Charge conservation 

At Q — the quark-photon vertex is completely spec- 
ified by the differential Ward identity 



iT b l {p,p;0) = ^-S^ 1 (p) 



_d_ 

dpti 



(13) 



If this is inserted in Eq. 
integration variables k — 



), one finds after a change of 



i 1 



Mt bb {P,Q)=2P ll F all {Q)=2N c 



A j4 



d 4 q 



Tr [f a V, q; P) S a (q) T a \q, q'; P) dS "^ p P) ] 



, (14) 



with q' = q — P. Comparing this expression with 
Eq. (^) with rj = 0, we recognize that the physical re- 
sult F(Q 2 = 0) = 1 follows directly from the canonical 
normalization condition for T ab with a BSE kernel K in- 
dependent of the meson momentum P. For the ladder 
truncation of the kernel, which we consider in our calcu- 
lation in the next section, this is the case. 

With a general momentum partitioning parameter r], 
the relation between the normalization condition and 
electromagnetic current conservation is not so obvious. 
However, using a different r/ in loop diagrams (without 
external quark lines) is equivalent to a shift in integration 
variables. For processes that are not anomalous, loop in- 
tegrals are independent of a shift of integration variables, 
provided that such a shift is performed consistently, and 
that all approximations employed respect Poincare in- 
variance. In performing such a shift, one has to take spe- 
cial care of the BSAs. The vertex function T(q, q'; P), as 
function of the incoming and outgoing quark momenta, 
does not depend on 77; it is only in commonly used decom- 
positions in terms of Lorentz invariant amplitudes such 
as Eq. (j8|), where 77 becomes relevant. The amplitudes E, 
F, G, and H are scalar functions of k 2 and k ■ P, which 
do depend on the choice for rj. Under a change of r], some 
of the different Dirac structures (e.g. the amplitudes F 
and G) will mix, as will the Chebyshev moments, fi in 
Eq. (M). Therefore, the results will be independent of 
the momentum routing in the loop integrals if and only 
if all Dirac amplitudes and their dependence on k ■ P 
are properly taken into account. Previously it has been 
shown that under these conditions the decay constants 
are indeed independent of 77 ||] . 

Use of a bare quark-photon vertex, in combination with 
dressed propagators, in Eq. (||), clearly violates charge 
conservation and leads to -^(0) 7^ 1. With the Ball-Chiu 
Ansatz [0, which is commonly used in DSE studies of 
electromagnetic interactions |lJ-0 ; the electromagnetic 
current is explicitly conserved, F(Q 2 = 0) = 1. However, 
the behavior of the form factor away from Q 2 = is not 
constrained by current conservation, and in the present 
model, use of the Ball-Chiu Ansatz leads to a value for r 2 
which is about 50% too small || . With the quark-photon 
vertex as the solution of the ladder BSE, together with 
quark propagators from the rainbow DSE, we satisfy all 
constraints from current conservation, and the calculated 
value of r 2 is within 5% of the experimental value ||. 
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D. Beyond rainbow-ladder truncation 

If one goes beyond the rainbow-ladder truncation for 
the DSEs for the propagators, BSAs and quark-photon 
vertex, one has to go beyond impulse approximation for 
the form factors in order to ensure current conservation. 
For example, one could include higher-order a s correc- 
tions to the rainbow-ladder DSE and BSE kernels, as 
depicted in Fig. JlL Following the general procedure de- 
veloped in Ref. p8[, one can show that both the WTI, 
Eq. (11), and the differential Ward identity, Eq. (|l3|), are 
preserved in the truncation indicated in Fig. |], as is the 
axial-vector WTI, which is important for the Goldstone 
nature of the pions. 





FIG. 1. The two leading-order vertex corrections to the 
rainbow DSE (top) and the corresponding five diagrams to 
be added to the ladder BSE kernel (bottom) for consistency 
with the relevant WTIs. The quark and gluon lines indicate 
dressed propagators in this and the subsequent figures. 

The resulting BSE kernel K(q,p; P) now becomes de- 
pendent on the meson momentum P, which means that 
the second term of the normalization condition, Eq.(g), 
is nonzero. To be specific, with the choice rj = 0, this 
introduces the four extra terms in the normalization con- 
dition, diagrammatically depicted in Fig. g. These four 
additional diagrams can be generated from the BSE ker- 
nel in the bottom part of Fig. [j] by taking the derivative 
with respect to the meson momentum P, where P flows 
through one quark propagator only. 







FIG. 2. The four diagrams due to the P-dependence of the 
kernel in the normalization condition, Eq. (g), if one includes 
the diagrams of Fig. ginto the DSE dynamics, and chooses all 
of the meson momentum P to flow through one quark only. 
The derivatives with respect to P are marked by slashes. 

Since taking the derivative with respect to P is equiva- 
lent to the insertion of a zero-momentum photon accord- 



ing to the differential WTI, Eq. (|11|), it is obvious which 
diagrams have to be added to the impulse approximation 
to ensure current conservation, see Fig. (g). In the limit 
Q — > these four additional diagrams become identical 
to the four additional diagrams in Fig. (g), provided that 
the vertex satisfies the differential WTI. Of course, there 
are similar contributions to A aba , which can be identified 
with terms in the normalization condition with r/ = 1. 







FIG. 3. The four additional contributions to the impulse 
approximation, Eq. (^), in order to ensure current conser- 
vation if one includes the diagrams of Fig. |l| into the DSE 
dynamics. 

Also, simple addition of contributions due to pion and 
kaon loops to Eq. (g), in combination with a ladder- 
rainbow truncation for the DSEs, will generally violate 
current conservation. Current conservation requires a 
consistent treatment of the kernels for both the DSE and 
BSE equations and the approximation for the photon- 
hadron coupling. At present it is not clear how to incor- 
porate meson loops self-consistently in such an approach, 
but we expect corrections coming from such loops to be 
small in the spacelike region. In Ref. (l9| it was demon- 
strated that the quark core can generate most of the pion 
charge radius, and that pion loops contribute less then 
15% to r 2 . For larger values of Q 2 the effect from me- 
son loops reduces even further, and for Q 2 > 1 GeV 2 we 
expect the contribution of such loops to be negligible. 



III. MODEL CALCULATIONS 



For the BSE we use a ladder truncation 



K(p,q;P) 



-g{k 2 )D^(k)^ 



7m ® ^Tlv : 



(15) 



where D^° c (k = p — q) is the free gluon propagator in 
Landau gauge. The resulting BSE is consistent with 
a rainbow truncation Y1(q,p) — > 7„A a /2 for the quark 
DSE, Eq. (g), in the sense that the combination produces 
vector and axial-vector vertices satisfying the respective 
WTIs. In the axial case, this ensures that in the chi- 
ral limit the ground state pseudoscalar mesons are the 
massless Goldstone bosons associated with chiral sym- 
metry breaking [g||. In the vector case, this ensures 
electromagnetic current conservation. 

The model is completely specified once a form is cho- 
sen for the "effective coupling" G(k 2 ). We employ the 
Ansatz g,[l(J 



4 



Qjk 2 

k 2 



+ 



47T 2 7m ^(fc 2 ) 



A. Results for uuu, usu, and uss form factors 



±ln 



1 



fc 2 /A 2 



QCD 



(16) 



with 7m = 12/(33 - 2N f ) and JF( S ) = (1 - exp ■ 4m . ; 
This Ansatz preserves the one-loop renormalization 
group behavior of QCD, and ensures that we reproduce 
perturbation theory in the perturbative region. The 
first term of Eq. ( |l6| ) implements the strong infrared en- 
hancement in the region < k 2 < 1 GeV 2 which is a phe- 
nomenological requirement for sufficient dynamical chiral 
symmetry breaking to produce an acceptable strength 
for the quark condensate (^0|. We use m t — 0.5 GeV, 
r = e 2 — 1, Nf = 4, Aqcd = 0.234 GeV, and a renormal- 
ization point (i = 19 GeV, well in the perturbative re- 
gion P,|l0||. The remaining parameters, us — 0.4 GeV and 
D = 0.93 GeV 2 , are fitted to give a good description of 



the chiral condensate, 



l>w/K 



and fn- The subsequent 



values for fx and the masses and decay constants of the 
vector mesons p, <p, K* are in agreement with the exper- 
imental data Jl(| , see Table |. 

TABLE I. Overview of the results of the model for the 
meson masses and decay constant, adapted from Refs. flp"o|l. 
The experimental value for the condensate is taken from 
Ref. [Eil. 





experiment p2[ 


calculated 




(estimates) 


C fitted) 


^u — cL 
"V=lGeV 


5 - 10 MeV 


5.5 MeV 


m /J=lGcV 


100 - 300 MeV 


125 MeV 


- mi 


(0.236 GeV) a 


(0.241 T ) a 




0.1385 GeV 


0.138 1 " 


U 


0.0924 GeV 


0.093 f 


rriK 


0.496 GeV 


0.497 1 " 


fK 


0.113 GeV 


0.109 


m p 


0.770 GeV 


0.742 


fp 


0.216 GeV 


0.207 


m K + 


0.892 GeV 


0.936 


Ik* 


0.225 GeV 


0.241 


mj, 


1.020 GeV 


1.072 


/* 


0.236 GeV 


0.259 



The pion and kaon form factors are given by 
FAQ 2 ) = lF udu {Q 2 ) + \F uM {Q 2 ) , 
F K +{Q 2 ) = fiWQ 2 ) + \F U - S - S {Q 2 ) , 
F K o(Q 2 ) = -\F dsd {Q 2 ) + \F dss {Q 2 ) , 



(17) 
(18) 
(19) 



where the quark and anti-quark charges are evident. We 
work in the SU(2) isospin limit, where the strong inter- 
action does not discriminate between u- and d-quarks, so 
for the pion we simply have Ft T {Q 2 ) — F uiiu (Q 2 ). Thus 
there are only three independent form factors, F uilu (Q 2 ) 1 
F u su(Q 2 ), and F uSS (Q 2 ), which are shown in Fig. [|. Our 
estimate of the numerical error in these calculations is 
less than 1% for F U u U (Q 2 ), and 2% for the other two 
form factors. 




0.05 0.1 

Q 2 [GeV 2 ] 

FIG. 4. The three independent form factors F U uu, F u su 
and F U ss- 

For the pion we use only the leading terms of the ex- 
pansion of the BSAs and the quark-photon vertex in k-P, 
Eq. (0). Higher order terms do not change the results 
more than 1% in this momentum regime, although they 
are needed at larger values of Q 2 . For the kaon we have to 
use more terms in the expansion, even at Q 2 = 0, to ob- 
tain independence from the parameter r), and to ensure 
current conservation. With terms up to order (k ■ P) 1 
only, there is a spread in our results of more than 10% at 
Q 2 = (from 0.94 to 1.06) if we change rj between and 
1. Including the next two terms reduces that spread to 
less then 3%, illustrating that the result of a loop integral 
is independent of this unphysical parameter 77, provided 
that all relevant Dirac structures and the dependence on 
k ■ P are properly taken into account. 

The results for F uau and F uSu are remarkably close 
to each other, indicating that the flavor of the specta- 
tor quark matters very little. Within our numerical er- 
rors, they are almost indistinguishable on the Q 2 domain 
shown. There is a slight difference in the slope of these 
form factors: r 2 au — 0.45 fm 2 versus r 2 Su = 0.47 fm 2 . 
These results are in good agreement with the pion charge 
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radius, r 2 = 0.46 fm 2 , obtained in Ref. |j] using all eight 
Dirac amplitudes of the quark-photon vertex. 

The result for F uSS is quite different in that it has a 
significantly smaller slope characterized by a radius pa- 
rameter r 2 ss = 0.21 fm 2 . This is due to the larger mass of 
the strange quark, and as a consequence the neutral kaon 
charge radius rL will be negative. A similar effect was 
observed for the neutron form factor, where the heav- 
ier mass of the + (wci)-diquark compared to the d quark 
mass leads to a negative charge radius |l7j . Our result is 
also consistent with the qualitative aspects of the vector 
meson dominance [VMD] picture: the lowest-mass bound 
state pole in the ssj- vertex is the <f>, at Q 2 = — 1 .0 GeV 2 , 
which is significantly further from the photon point than 
is the p pole in the uwy- vertex at Q 2 = —0.6 GeV 2 . This 
observation, as well as the difference between r 2 ss and 
r usm i s consistent with the larger mass of the strange 
quark. 



Finally, in Fig. ^ we present Q 2 F{Q 2 ) for ir and K 0,± 
for a larger Q 2 range to anticipate data that may be 
forthcoming from experiments at JLab ||Q and possibly 
other facilities in the future. In this momentum range, 
even for F 1T {Q 2 ) the dependence on k ■ P becomes impor- 
tant, and terms up to (k ■ P) 3 in Eq. @ are required to 
produce a converged result at Q 2 = 1 ~ 3 GeV 2 . Higher- 
order terms do not change the results by more than 1% 
in this momentum range. Our estimate is that the net 
numerical accuracy for F uau , F uSu , and F uSS is about 2- 
3% at these values of Q 2 . This translates to a similar 
level of accuracy for F n and F K + , and to a somewhat 
larger relative error, about 5%, for F K o, which is the dif- 
ference of F uSS and F uSu . At Q 2 > 3 GeV 2 , higher-order 
Chebyshev moments may be necessary, but current nu- 
merical methods prevent their accurate determination at 
large Q 2 . 



B. Results for the meson form factors 

The results in this model for the pion form factor at low 
Q 2 , in particular the pion charge radius, were presented 
previously |^| . The obtained charge radii for the kaon are 
presented in Table |l[ and are in reasonable agreement 
with the experimental data, without any readjustment of 
the model. In Fig. || we show our result for the charged 
kaon, which is in good agreement with the available data. 



1.0 



0.9 



O 0.8 



0.7 



0.6 




0.04 



0.08 0.12 
Q z [GeV 2 ] 



0.16 



FIG. 5. The calculated K + form factor compared to the 
data from Refs. J26) (open diamonds) and (solid squares). 
Within numerical errors, |-Fa'+(0)| 2 = 1- 



TABLE II. Our results for the charge radii, compared with 
the experimental values given in Refs. [E3l— E5l| . 



charge radii 


experiment 


calculated 


r 2 


0.44 ± 0.01 fm 2 


0.45 fm 2 


r 2 


0.34 ± 0.05 fm 2 


0.38 fm 2 


r K° 


-0.054 ± 0.026 fm 2 


-0.086 fm 2 




K + , up to 1 51 order in Eq. (9) 

K + , up to 3 rd order 

— K° 




2 2.5 
Q z [GeV 2 ] 

FIG. 6. Q 2 times the kaon form factors (top) and pion 
form factor (bottom). The pion data are from Ref. p7| . 

Over the entire spacelike momentum range considered, 
Fn(Q 2 ) < F K + (Q 2 ), and Q 2 F(Q 2 ) rises with Q 2 until 
Q 2 = 3 GeV 2 for all three form factors. In this mo- 
mentum range our results for both the pion and the K + 
form factor can be fitted quite well by a simple monopole 
m 2 /(Q 2 + m 2 ), with a mass m 2 — 0.53 GeV 2 for the 
pion and m 2 = 0.61 GeV 2 for the K+ . A VMD model, 
two monopoles for the two form factors F uSu and F uSS 
in Eqs. (18) and ( |l9| ) with the physical p and <f> masses 
respectively, does not reproduce our results for the kaon 
form factors very well. For example, at Q 2 = lGeV 2 , 
VMD overshoots our F K + calculation by almost 10%, 
whereas the monopole fit is within 2% of our result. This 
difference between VMD and our calculations grows with 



6 



Above Q 2 ~ 3.5 GeV the monopole fits begin to de- 
viate significantly from our results and Q 2 F(Q 2 ) starts 
to decrease. In a more realistic model, that takes meson 
loop corrections into account self-consistently, it is very 
well conceivable that this turn-over happens at somewhat 
lower values of Q 2 : meson loops are expected to con- 
tribute up to 15% to r 2 fll9f , but their contribution to 
the form factor decreases rapidly with increasing space- 
like momenta. In the presence of meson loop corrections 
the contribution to the form factor from the impulse ap- 
proximation has to be smaller than in our calculation 
in order to maintain agreement with the low-Q 2 data. 
Therefore it is not unlikely that at intermediate momenta 
in the present approach we overestimate the form factors, 
which may explain the difference between the data points 
at Q 2 = 3.3 GeV 2 and our calculated results. More ac- 
curate results from JLab, in combination with realistic 
model calculations that include meson loop corrections 
self-consistently, may be able to resolve this question. 

At asymptotically large Q 2 , factorized pQCD Kg] pre- 
dicts that the form factor behaves like Q 2 F(C^) — > c, 
with c— 16-7T f 2 a s {Q 2 ). Since our truncation and the 
Ansatz, Eq. (16), is constructed so as to preserve asymp- 



totic freedom, we are guaranteed to recover the leading 
power-law asymptotic behavior. An explicit verification 
of this behavior, and calculation of the constant c, is not 
readily available within our present framework since nu- 
merical accuracy at large Q 2 is problematic. However, it 
is clear from our results that at Q 2 ~ 4 GeV 2 the form 
factor has not yet reached its asymptotic value, and it is 
unlikely that experiments can access the true asymptotic 
region in the near future. In simplified models such as 
that of Ref. jL5j however, it is straightforward to demon- 
strate that the impulse approximation does indeed lead 
to the power-law behavior predicted by pQCD. 



the available data over the entire Q 2 range considered, 
and the calculated charge radii are within the error bars 
of their experimental values. 

These charge radii are somewhat larger than those ob- 
tained in a previous study [l3| that was framed in terms 
of semi-phenomenological representations for BSAs and 
confined quark propagators within the impulse approxi- 
mation. The main difference with that work is that here 
we use numerical solutions of truncated DSEs for all the 
elements needed in Eq. (Q), and that all our parame- 
ters were fixed previously. In comparison with theoreti- 
cal calculations based on other methods, it is interesting 
to note that our results are very similar to those obtained 
in Ref. |2^], in particular for the neutral kaon. 

At intermediate values of Q 2 our calculations are qual- 
itatively similar to those obtained in both Ref. JrJ 
and Ref. §§]. Up to about Q 2 = 3GeV 2 , both F v 
and F K + can be fitted quite well by a monopole 
form, with monopole masses of m 2 = 0.53 GeV 2 and 
to 2 = 0.61 GeV 2 respectively. At large Q 2 the DSE 
approach does reproduce the pQCD power-law behav- 
ior pg|, but this behavior does not occur until well be- 
yond! 01 the Q 2 range considered in our present calcu- 
lations and accessible at current accelerators. 
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IV. SUMMARY 



We calculate the pion and kaon electromagnetic form 
factors within the DSE approach. The method is com- 
pletely Poincare invariant, and the only approximation 
made is a self-consistent truncation of the set of DSEs, 
which respects the relevant vector and axial- vector WTIs. 
The employed quark propagators, the meson BSAs, and 
the quark-photon vertex are solutions of their DSEs in 
rainbow-ladder truncation with all parameters fixed pre- 
viously by fitting the chiral condensate, m^/K and / w . 
We include all relevant Dirac amplitudes for the BSAs 
and their dependence upon k ■ P. The electromagnetic 
current is explicitly conserved in this approach, and there 
is no fine-tuning needed to obtain F 7T (Q) = 1 = F K +(0) 
and Fjc ° (0) = 0. We also demonstrate explicitly that our 
results are (within numerical accuracy) independent of 
the momentum partitioning of the BSAs. The obtained 
pion and kaon form factors are in good agreement with 
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